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ABSTRACT 

We present the results of a high-resolution X-ray imaging study of the stellar population in the 
Galactic massive star- forming region RCW49 and its central OB association Westerlund 2. We 
obtained a ~40 ks X-ray image of a ^17'xl7' field using the Chandra X-ray Observatory and deep 
near-infrared (NIR) images using the Infrared Survey Facility in a concentric '--^S.'SxSfS region. We 
detected 468 X-ray sources and identified optical, NIR, and Spitzer Space Telescope mid-infrared 
(MIR) counterparts for 379 of them. The unprecedented spatial resolution and sensitivity of the 
X-ray image, enhanced by optical and infrared imaging data, yielded the following results: (1) The 
central OB association Westerlund 2 is resolved for the first time in the X-ray band. X-ray emission is 
detected from all spectroscopically-identified early-type stars in this region. (2) Most (~86%) X-ray 
sources with optical or infrared identifications are cluster members in comparison with a control field 
in the Galactic Plane. (3) A loose constraint (2-5 kpc) for the distance to RCW49 is derived from 
the mean X-ray luminosity of T Tauri stars. (4) The cluster X-ray population consists of low-mass 
pre-main-sequence and early-type stars as obtained from X-ray and NIR photometry. About 30 new 
OB star candidates are identified. (5) We estimate a cluster radius of 6'-7' based on the X-ray surface 
number density profiles. (6) A large fraction (~90%) of cluster members are identified individually 
using complimentary X-ray and MIR excess emission. (7) The brightest five X-ray sources, two Wolf- 
Rayet stars and three O stars, have hard thermal spectra. 

Subject headings: X-rays: stars — infrared: stars — ISM: H II regions: individual (RCW49) — open 
clusters and associations: individual (Westerlund 2) — stars: Wolf-Rayet — stars: 
pre-main-sequence 



1. INTRODUCTION 

H II regions, which are recognized by diffuse hydrogen 
emission ionized by UV photons from stars earlier than 
B2, are the sites of massive star formation. Strong ra- 
diation and winds from early-type stars in massive star- 
forming regions (SFRs) have effects that cannot be un- 
derstood by a simple extrapolation of low-mass SFRs. 
X-ray studies of massive SFRs, which are potentially im- 
portant to detect hot gas as well as stellar constituents, 
had lagged behind those of low-mass regions in the past 
due to their greater distance, source confusion, and ob- 
scuration. Both high spatial resolution and hard X-ray 
sensitivity are required to overcome these difficulties. 

The Chandra X-ray Observatory is capable of conduct- 
ing high-resolution and high-sensitivity imaging observa- 
tions in the 0.5-8.0 keV band. Its ~1" spatial resolution 
and sensitivity are sufficient to resolve scales ^10^^ cm 



(a typical size of ultra-compact H II regions; Church 



well 



2002 1 and to detect X-ray emission from embedded 



late-type pre-main-sequence sources at a distance of sev- 
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eral kiloparsecs, typical of well-studied Galactic massive 
SFRs. 



Along with the Chandra Orion Ultradeep Project ( Get- 



man et al. 2005 1 , which led to a series of X-ray imag- 
ing studies of the nearest {D < 0.5 kpc) rich SFR, we 
have seen a surge of results from X-ray studies of mas- 
sive SFRs at larger distances. Galactic massive SFRs 
with X-ray imagin g surveys of their stellar contents in- 
clude : NGC3603 ( jMoffat eia l.120021), W3 (IHofner et al 



20021), the Carina ne bula (|Evans et al. 2003, 2004; San- 



chawala et al.||2007|), NGC^30 (jPamiani et al. 2004 
the T rifid Nebula (|Rho et al.|2004|^ 106 ( |Giardino et al 
2004t, Ara OBI ( [Skinner et_a r 2005iCe£heus B ( |Ue¥ 
mail et al.|2006[ ), RCW 3 8 ([Wolk et all2006l) , the ArcEes 
and Quintuplet clus ters dWang et al.pOOGi, NGC 2362 
d Pelgado et al.|2006|, nGC'l 'lU fVlaccomio et al.'2006j, 
Westerlund 1 ( |skinn er et al. 2006), IC 1396N ( Getmanl 
et al.|2007aL NUIT5S57 (Wang et al. 2007al, the Rosette 



et al. 2007a1), NUIT5 '357 (Wang et al. 2007a |), the R osette 
nebula ([Townsle^t al.,,2003, Wang e t al.'^2007b') , M 16 
(Lins ky^et al.||2007D , Cygnus 0B2 ([Al bacete Colo mbo et 
al. 2007]), and M17 ([Townsley et aT 2003 Broos eTaT 



2007p, andlvm ( [Townsley 
i 20C)7j). Early result s are reviewed in Townsley (2006b] 
and jFeigelson et al"] (]2007[). 



This paper describes the results of the Chandra imag- 
ing study of the stellar contents in a southern massive 
SFR RCW 49 and its central OB association Westerlund 
2. The region was imaged in the mid-infrared (MIR) 
band as a part of the Galactic Legacy Infrar ed Mid-Plane 
Survey Extraordinaire program (GLIMPSE; 'Benj amin et 
al. 2003) using the Spitzer Space Telescope (Werner et 
al. 2004). We also obtained deep near-infrared (NIR) 
images of this region using the Infrared Survey Facility. 
Many previous Chandra studies were combined with NIR 
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imaging data, but none was presented with MIR imag- 
ing data, which is a growing tool for the stellar census 
studies in massive SFRs with the launch of Spitzer. We 
demonstrate the abihty of high-resolution X-ray imaging 
to identify and understand the nature of the stellar pop- 
ulation in massive SFRs with a particular emphasis on 
the combination with NIR and MIR imaging. 

The outline of this paper is as follows. In § |2] we 
give a brief review of RCW49 and Wester lund 2. In 
§ |3] we present our X-ray and NIR observations of this 
region. In § [4] we describe our X-ray and NIR data 
reduction and the identification of X-ray sources with 
optical, NIR, and MIR sources using our deep NIR data 
as well as archived data. In the data reduction, special 
care is taken in source extraction, photometry, and cross- 
correlation in crowded regions. In total 468 X-ray sources 
are detected and 379 of them are identified by optical and 
infrared sources. 

In § [5] we discuss how the X-ray data are relevant 
to discriminate intrinsic RCW49 members, identify new 
sources, and infer their nature. Seven multi-faceted and 
loosely-related results ar e dis cussed in each subsection 
In the first subsection ( 
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and 



Moffat et al. ( 1991 ) identified six 06-7 stars based 



on low -resolution optical spectroscopy. Rauw et al 



(20071 added six more 0-type stars and gave a con- 
sistent classification based on medium-resolution optical 
spectroscopy. These stars were reclassified towa rd earlier 
spectral classes ranging between 03 and 06.5 (Rauw et 
p. 2007) . 

I'hree more early-type stars are found beyond the clus- 
ter core. O ne is a star with a spectral type of 04V(f) 



or 05V(f) dUzpen et al][2005| |Rauw et al.| 2007|, which 
we refer to as MSP 18 (IjVlc 



ottat et al.lll99li. Another is 



a binary Wolf-Rayet star WR 2 0aTWN6-HWN6; van der 



Hucht 2001 Rauw et al. 2005 



5.1 1, we show that the central sive binary in~the Galaxy wit. 



OB association is resolved in the X-ray image and all the 
massive members in this region are detected in the X-ray. 
In § |5.2| we argue that most of the X-ray sources with 
optical or infrared counterparts are intrinsic to RCW 49 
by a comparison study with the X-ray population of a 
control field in the Galactic Plane. In S 15.31 we obtain 



which is the most mas 
a well-determined mass 
Rauw et al.||2005|. Another Wolf-Rayet s tar WR20b 
van der Hucht||2001|) lies ~3f7 from Wd2 CShara et al 



a loose constraint on the distance to RCW 49 using the 
mean X-ray luminosity of T Tauri stars in the mass range 
of 2.0-2.7 Mq. 



In 



5.4 we evaluate the nature of X-ray 
sources using X-ray and infrared photometric data. We 
employ the X-ray versus NIR flux and luminosity plots 
in addition to the conventional NIR color-magnitude di- 
agram. We argue that the majority of X-ray sources 
are low-mass pre-main-sequence sources and early-type 
main sequence stars. Some early-type stars are found 
to be harder and more luminous in the X-ray than oth- 
ers with a similar NIR brightness. In § |5.5| we study 
the spatial distribution of X-ray-identified cluster mem- 
bers. Because of the low level of contamination of X-ray 
samples, we can identify cluster members far from the 
central OB association. We find a possible excess pop- 
ulation at ~4' and the cluster boundary at 6'-7' from 
the cluster center. In § 5.6 we examine the fraction 



of RCW 49 members individually selected using X-ray 
emission and find that about half of the RCW 49 sources 
brighter than K^^U mag are identified by X-rays. The 
rate increases to ~90% when combined with MIR excess 
emission, another indicator of cluster membership com- 
plementary to the X-ray emission signature. In the last 
subsection (§ 5.7 1, we discuss the X-ray spectra of the 
five brightest X-ray sources. We summarize these results 

2. RCW 49 AND WESTERLUND 2 

RCW 49 was identi f ied as an optically-visible H II re- 
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2004 Rauw et al. 2007 Ascenso et 



lar constituents m t he region were 



studied by imaging studies in the optical ( Moffat & Vogt 



991[ ). I'hese early-type stars in and around Wd2 are 
probably the cause of exotic high energy phenomena dis- 
cover ed in this region^ such as extended hard X-ray emis- 
sion (Townsley et al. 20041 found in the same Chandra 
data set presented in th is paper and TeV gamma-ray 
emission (Aharonian et al. 2007) found by the High En- 
ergy Stereoscopic System telescope. 

A much larger number of lower-mass members are an- 
ticipated, many of which are at the pre-main-sequence 
stage that can be identified either by infrared excess or 
X-ray emission, arising respectively from circumstellar 
disks and envelopes and from enhanced coronal magnetic 
activity. Some of the low-mass members were detected 
by mid-infrared exce ss emission using the Spitzer data 
( Whitney et al.|2004 V In the X-ray band, the two previ- 
ous studies did not have sufficient spatial resolution and 
sensitivity to distinguish individual stars. In these stud- 
i es, respectively using Einstein a nd ROSAT, only one 
(|G"oldwurm et al.|1987 l and seven (Belloni & Mereghetti 
1994[ ) X-ray sources were detected in the Ctiandra stud- 
ied held presented here. Most of these sources are iso- 
lated early-typ e stars and unresol ved stars in Wd2 (IE 
1022. 1-5730 inlGold wurm et al.'1987 and source number 
24 in [Belloni fc Mer eghetti 1994). 

The ionized gas content at ~10'* K in RCW 49 was 
studied by a radio con tinuum imaging study at 0.843, 
1.38, and 2.38 GHz ( Whiteoak fc Uchida|[l997t . Two 
shells appear in the images, one around the complex of 
Wd 2 and WR20a with a diameter of ~7f3 (radio ring A) 
and the other around WR 20b with a diameter of ~4fl 
(radio ring B). Enhanced emission was found at the po- 
sition where these two shells overlap. A similar global 
structure was found in the dust content. The MIR im- 
age taken by the Infrared Array Camera (IRAC) onboard 
Spitzer traces the dust distribution via polycyclic aro- 
matic hydrocarbon emis sion at its three ([3.6], [5 .8], and 
[8.0] ^m) of four bands ( |Churchwell e t al. 2004). More- 
over, the Spitzer image reveals a network of highly struc- 
tured filaments, pillars, and shocks at arcsecond scales. 
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some of which are suggestive of instabihties. 

Both gas and dust distributions are dominated by the 
massive stars in Wd 2 and the two Wolf-Rayet stars via 
their strong stellar winds and UV radiation that sculpt 
the surrounding interstellar medium. These effects not 
only suppress further star formation, but also trigger the 
second generation of star formati on. From the s patial 
distribution of MIR excess sources, Whitney et al. ( 2004 1 
argued that star formation is taking place preferentially 
in ^i' around Wd 2. The radius corresponds to the outer 
boundary of the radio ring A where swept-up material is 
accumulated, suggesting that these young sources are a 
consequence of induced star formation. 

The distance to RCW 49 is highly c ontroversial, rang- 
ing fr om 2 to 8 kpc in the literature (Churchwell et al. 



2004] and references therein). The kinematic distance de- 
termination gives a poor constraint because the region is 
in the tangential direction of the Saggitarius-Carina arm. 
The de bate on the distanc e continues in two most recent 
papers; Rauw et al. (20071 argued for 8.0±1.4 kpc based 
on optical s pectro-p notometric study of bright stars in 
Wd2, while Ascenso et al. (2007) proposed 2.8 kpc based 



on NIR magnitudes and colors of RCW 49 sources on the 
Henyey track. With the present X-ray data, we can also 
loosely constrain the distance using the mass-stratified 
mean X-ray luminosity (§ |5.3| . This gives an estimate 
of 2-5 kpc. Despite intensive study, the distance to this 
cluster remains an open questio n. We adopt 4.2 kpc as 



used by [Churchwell et al.| ( |2004[ ) , which is roughly a ge- 
ometric mean of the smallest and the largest distance 
estimates, to avoid the worst case when the question is 
settled. 

The association of Wolf-Rayet stars and an 04-05 V(f) 
star with Wd2 suggests that the star cluster is no 
older than a few Myrs. No cluster member appears 
to be located away from the zero-age-main-sequence 
isochrone c urve in the o ptical color-magnitude diagram, 
from which Carraro fc Munari (2004) estimated an age 
of < 2 Myr. Ascenso "et al., ( ^2007j similarly derived the 
cluster age of ~2 Myr by fitting the color-magnitude dis- 
tribution of NIR sources with theoretical calculations. 

3. OBSERVATIONS 

We carried out an X-ray observation of RCW 49 us- 
ing the Advanced C CD Imaging Spectrometer (ACIS; 
Garmire et al. 2003 ) onboa rd the Chandra X-ray Ob- 



servatory ( ,Weisskopf et ar|l2002^ from 2003 August 23 
UT 18:20 to August 24 U 1 4:54. Four imaging array 
(ACIS-I) chips covered a 17'xl7' field centered at (R. A., 
decl.)=(10'^24'"00!5, -57°45'18") in the equinox J2000.0 
for a 36.7 ks exposure. ACIS-I covers the 0.5-8.0 keV en- 
ergy band with a spectral resolution of ~150 eV at 6 keV 
and a point spread function (PSF) radius of ~0'.'5 within 
~2' of the on-axis position, degrading to ~6" at a 10' 
off-axis angle. The data were taken with the very faint 
telemetry mode and the timed exposure CCD operation 
with a frame time of 3.2 s. 

We conducted near-infrared (NIR) observations on 
2004 December 25 and 28 using the Simultaneous three- 
color InfraRed Image r for Unbiased Surveys (SIRIUS; 
Nagayama et al.||2003 1 mounted on the Cassegrain focus 
of the Infrared Survey Facility (IRSF) 1.4 m telescope at 
the South African Astronomical Observatory. SIRIUS is 



a NIR imager capable of obtaining simultaneous images 
in the J, H, and Kg bands using two dichroic mirrors. 
The instrument is equipped with three HAWAII arrays 
of 1024x1024 pixels. The pixel scale of 0'.'45 is an excel- 
lent match with the on-axis spatial resolution of Chan- 
dra. With ten-point dithering, we covered 8f3x8f3 fields 
at two positions, one aimed at RCW 49 (10^24™01!9, - 
57°45'31") and the other at a control region (10'^27'"13!5, 
-58° 01 '26") free from very bright sources displaced by 
~30' from RCW 49. Both fields have the same Galactic 
latitude. 

Each frame was exposed for 30 s and the total expo- 
sure time was ~30 minutes equally for the two regions. 
Additional short (5 s) exposure frames were also taken 
for the two regions to extend the dynamic range of the 
data set to brighter sources. The observing conditions 
were photometric with a seeing of < I'.'O on both nights. 

Figures fl] (a) and (b) respectively show the ACIS and 
SIRIUS images of the study field. The SIRIUS image has 
a wid er coverage by ~2. 9 times than the image presented 



in [Ascenso et al. (20071 



4. DATA REDUCTION AND ANALYSIS 

4.1. X-ray Data 

4.1.1. Reduction & Source Extraction 

We reprocessed the raw data distributed by the Chan- 
dra X-ray Center to obtain an X-ray event list. A back- 
ground rejection algorithm specific to data taken with 
the very faint mode was applied. Events were further 
cleaned by removing cosmic-ray afterglow and applying 
filters bas ed on the event grades , status, and good-time 



intervals (Townsley et al. 2003|). Charge transfer inef- 



ficiency was corrected fltawnsley et al. 2002 1 and the 



positions of all events were improved usmg a sub-pixel 



repositioning technique (Tsunemi et al. [2001 1 

We detected sources using the wavdetect algorithm 
in the CIAO package* independently for the soft (0.5- 
2.0 keV), hard (2.0-8.0 keV), and full (0.5-8.0 keV) band 
images. The software failed to detect sources in crowded 
regions, where we manually added candidate sources with 
the aid of smoothed images in different kernel sizes and 
PSF-deconvolved images obtained using a maximum like- 



lihood technique ( Lucy[[l974 Townsley et al 



2006a I. In 



total, 556 X-ray source candidates were identihed includ- 
ing 60 manually added sources. We took a stance to add 
sources manually in a liberal manner here, and we exam- 
ined their validity in a rigorous manner later. 

For all the candidate sources, we used ACIS Extract^ 
version 3.107 for automated systematic source and back- 
ground event extractions. The source positions were 
refined by correlating the event distribution with the 
PSF or by using the centroid of the events. Signals 
were then extracted from a 90% encircled energy poly- 
gon of 1.5 keV X-rays to derive source counts in 0.5- 
8.0 keV (Csic)- When two extraction regions overlap, 
the regions were reduced. A masked dataset was cre- 
ated by removing the events in a circle with a ^--^99% 
encircled energy radius around all sources, from which 



Seejhttp://asc. harvard .edu /ciao/ 1 for detail. 
See the ACJlb Extract User's 



Guide at 



[http:// www.astro.psu.edu/xray/docs/TARA/ae_users_guide Thtinll 
lor detail. A complete descriptio n ot the procedure usi ng the A(Jis 
Extract package can be found in [Getman et al.| ( [2005^ . 
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Fig. 1. — (a) Smoothed X-ray image. Red, green, and blue indicate the intensity in the 0.5-1.7, 1.7-3.8, and 3.8—8.0 keV bands, 
respectively. The fields of view of the four ACTS chips and SIRIUS ar e shown by solid and dash ed squares. Free-style lines taken from 
[Churchwell et ah] Ij2004j) indicate the global features in the radio image ( [Whiteoak fc Uchida|1997l l. (b) NIR image. Red, green, and blue 
indicate the mterisity in the Ks, H, and J bands, respectively. The positions ot some landmark sources are given in small squares in both 
images and labeled in panel (a). 



background events were extracted locally around each 
source to derive their background and net counts, Cbkg 
and Cnet = Csrc " Chkg X y4srcMbkg, where As^c and 
^bkg are the integrals of the masked exposure map val- 
ues within the source and background extraction regions. 

Each source candidate was tested for validity using 
the photometric significance (PS) and probability of no 
source (Pb) statistics extracted by ACIS Extract. PS 
is the commonly used photometric signal-to-noise ratio 
given by 

PS=^, (1) 



where the approximation by Gehrels (1986) was em- 
ployed to calculate the uncertainty AC'net- -Pb is the 
probability for the null hypothesis that all detected 
counts or more are explained by fluctuations of the back- 
ground and is defined by 



Pb = 1 



E 

n=0 



P(Cbkg X ^,cMbkg,n), (2) 



where P(A,n) indicates a Poisson probability distribu- 
tion function of a mean A to have n counts. 

We recognize 468 source candidates with PS > 1.0 and 
Pb < 1.0x10"^ to be valid X-ray sources (Table]!]). None 
of them is bright enough to cause pile-up. Almost all the 
sources are newly detected here. 

4.1.2. Photometry & Time Vanability 

For all valid sources, ACIS Extract was used for re- 
calculating their source and background counts, produc- 
ing instrumental responses, and constructing spectra and 
light curves for systematic photometry and spectroscopy 
analyses. The photometry results are compiled in Ta- 
ble m Columns (1) and (2) are the sequential number 
and the source name. Columns (3)-(6) give position 
information (R. A. and decl. in the equinox J2000.0, 
their uncertainty, and the off-axis angle). Columns 



(7)-(9) show counts in 0.5-8.0 keV; Cnet, ACnct, and 
^bkg — ^bkg X AsrcMbkg, whereas column (10) shows 
the net counts in 2.0-8.0 keV (Cnot.hard)- Column (11) 
indicates the PSF fraction of source extraction, in which 
smaller numbers than the default value (~0.9) indicate 
that the extraction region was reduced to avoid overlap- 
ping. Columns (12) and (13) are for PS and Pb. Columns 
(14) and (15) are sets of flags to indicate extraction 
anomalies and source variability. Column (16) gives the 
effective exposure time, the time needed to accumulate 
the observed counts if the source were at the aim point. 
Column (17) lists the median energy of the extracted 
X-ray events corrected for the background. Median en- 
ergy is proposed to be a better quantity to characterize 
X-ray spectral hardness and to give a reliable estimate 
of extinction than the conventionally-us ed hardness ratio 
( |Hong et al.|2004[|Feigelson et al.||2005|. Column (18) is 
the photometric flux discussed in jj |4.1.3| 

X-ray variability was examined for 377 bright sources 
off the chip gaps using a Kolmogorov-Smirnov test. The 
null hypothesis probability (Pks) that the flux is con- 
stant was tested and sources were placed into three 
classes; (a) no evidence for variability (Pks > 5 x 10~^), 

(b) possibly variable (5 x 10~^ < Pks < 5 x 10^^), and 

(c) deflnitely variable (Pks < 5 x 10""^). Seven sources 
are in the "deflnitely variable" class, showing flare-like 
light curves. 

4.1.3. Spectroscopy & Flux Estimates of Famt Sources 

Table 2] gives the results from the X-ray spectral flts 
obtained with ACIS Extract, where 230 spectra with 
PS > 2 were fltt ed by an optically -thin thermal plasma 
(APEC) model ( [Smith et aL]|2001| ) with interstellar ab- 
sorption. The best-flt values of plasma temperature 
{k^T), the volume emission measure (EM), and the ex- 
tinction column density (iVn) were derived. The pho- 
tometric and spectroscopic tables (Tables [l] and [2| are 
in a format similar to those found in associated papers 
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most of the sources in Table [2] For the remaining sources, 
which are mostly faint, we estimated their X-ray flux 
using broad-band photometry: the mean energy times 
the net count rate di vided by the mean effective area 
( Tsujimoto et al. 2005 1 . The consistency between the two 
flux estimates was checked for the brightest 173 sources 
(Fig. [2|, for which both values are available. They are 
overaUin good agreement with a systematic offset (20%). 
We hereafter use the spectroscopic flux for the 173 bright 
sources and the photometric flux corrected for the offset 
for the remaining sources. 
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Fig. 2. — Consistency between the two X-ray flux estimates in the 
0.5—8.0 keV band. The values derived from spectral fits (abscissa) 
are plotted against those from broad-band photometry (ordinate) 
for 173 sources with successful thermal fits in Table[2] The dashed- 
and-dotted line indicates y = x, while the dashed line indicates the 
best-fit model of y = 1.20a;. 



4.2. NIR Data 

All SIRIUS frames were reduced following standard 
procedures using IRAF, including dark-current correc- 
tion, sky subtraction, and flat-flelding. Dithered images 
were registered to construct combined images in the three 
bands (J, H, and Kg). Sources were extracted and their 
photometry was measured using the daoflnd and daophot 
tasks in IRAF. 

In order to cope with source crowding in the SIRIU S 
images, we used the following procedure ( Stetson||l987 l. 
First, we constructed the average PSF from dozens of 
bright sources across each image. The proflle of each ex- 
tracted source was fltted with the average PSF and the 
source was masked from the image if the flt was suc- 
cessful. The procedure was repeated until most sources 
were masked. We then reconstructed the PSF from the 
masked image, which we consider is freer from contami- 
nation than the original one. By repeating the procedure, 
we separated closely spaced sources and deblended their 
overlapping signals. 

We extracted 10,540 and 9,768 sources in RCW 49 and 



the control flelds, respectively, at levels larger than 3 
tr of the sky background noise in all three bands. We 
derived their positions and magnitudes in the J, H, and 
Ks bands in the longer-exposure images. We replaced 
magnitudes with those derived from the shorter-exposure 
images when they were saturated at levels brighter than 
11.0 (J), 11.5 (iJ), and 12.0 (Ks) mag. The 3 a detection 
limits of the longer-exposure images in the RCW 49 field 
are ~20.5 (J), ~19.8 (H), and ~18.4 (K^) mag, whereas 
the 10 cr limits are ~19.0 (/), ~17.8 (H), and -16.8 (/Q) 
mag. These limits were ~0.5 mag shallower than those 
measured in the control field because of the larger source 
confusion. All the detected sources are used for the X- 
ray counterpart search, while only those with magnitude 
uncertainty less than 0.1 mag are used for photometric 
analysis. 

4.3. Optical and Infrared Identifications of X-ray 
Sources 

We identified the X-ray sources in the optical, NIR, 
and MIR bands using the Naval Observatory Merged 
Astrometric Dataset (NOMAD) the Tw o-Micron All- 
Sky Survey (2MASS; 'Skruts kie et al.p97| , SIRIUS, and 
GLIMPSE. GLIMPSE is a Legacy Prograin of the Spitzer 
Space Telescope to survey a ~220 deg^ region of the 
Galactic Plane at four MIR bands ([3.6], [4.5], [ 5.8], and 
[8. 0] iim) with a 1'.'6-1'.'9 resolution using IRAC ([Fazio et 
ar]|2004[ ). A l?7x0?7 region encompassing RCW 49 was 
mapped with ten 1. 2 s exposures in the initial phase of 
the program (Church well et al. 2004). 

First, the ACIS image was shifted to match the ab- 
solute astrometry of 2MASS using the closest ACIS- 
2MASS pairs. The residual displacements of the pairs 
in the R. A. and decl. directions are plotted in Fig- 
ure ^ in which the rms (1 cr) of the displacements is 
~0'.'56 (Fig. [3] the inner dotted circle). We recognized 
216 ACIS-2MASS pairs within 2 a displacement (Fig. fs] 
the outer dotted circle) to be physical counterparts. We 
similarly found NOMAD, SIRIUS, and IRAC counter- 
parts within 2 cr of ACIS sources after a correction of 
the systematic offset, where 1 cr is ~0'.'63, ~0'.'38, and 
~0'.'59, respectively. Consequently, 230, 299, and 156 
ACIS sources were found to have NOMAD, SIRIUS, and 
IRAC counterparts. Both the numbers of false positives 
(unrelated pairs identified as counterparts) and false neg- 
atives (physical pairs unidentified as counterparts) are 
estimated to be <10. 

The resultant identifications and the optical (R), NIR 
(J, H, and Ks), and MIR ([3.8], [4.5], [5.8], and [8.0]) 
photometry are summarized in Table |3] The i?-band 
magnitude is from NOMAD, which is a compilation work 
of several catalogues. NIR magnitudes are either from 
2MASS or SIRIUS, which is indicated by the NIR flag 
column. 

In total, 379 of the 468 ACIS sources have either NO- 
MAD, 2MASS, SIRIUS, or IRAC counterparts. Except 
for massive stars in and around Wd2, almost all the 
X-ray sources are poorly studied in other wavelengths. 
Spectroscopic cla ssifications are available for five source s 
outside of Wd 2 (Moffat et al j 
HD 90273 (0 7), RD'iffmT 
MSP 158 (com), and MSP218 (AO III), among which 
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Fig. 3. — Astrometric consistency between ACIS and 2MASS 
after correcting for a systematic offset. The differences of R. A. 
and decl. between the ACfS and 2MASS positions (AR. A. and 
ADecl.) are plotted for all the closest pairs. The dotted circles are 
drawn at f and 2 a radius, whereas the solid square represents the 
ACIS pixel size of 0'.'492x0y492. 

only HD 90273 was detected in X-rays. The identifica- 
tions with sources in earher work are also given in Ta- 
ble El 

5. DISCUSSION 

5.1. Early-type stars in Westerlund 2 and Beyond 

Figure |4] shows a close-up view of Wd 2 in the SIRIUS 
and ACIS images. The Chandra observation of RCW49 
resolved X-ray sources in the central OB association. A 
dozen early-type stars from 07 to 03 are present in the 
cluster core ; MSP 151, 157, 167, 171 175, 182, 183, 188, 
199, 203, and 263 (Rauw et alT]|2007[ ). In addition, three 
early-type stars, an 04 V(frorT?5 V(f) star (MSP 18; 
Uzpen et al.|2i 
IWa 



2005 



Rauw et al.|2007| and two Wolf-Rayet 



stars (WR20a = MSP 240 and WR20b), he outside the 
cluster core. 

We detected X-ray emission from all of these early- 
type sources. The identifications are given in Table [3] 
The results illustrate clearly that high-resolution X-ray 
imaging observations, even with a short exposure of 30- 
40 ks, are powerful enough to resolve and identify the 
crowded massive members in OB associations at a few 



kiloparsecs aw ay, as was found by Wang et al. (2007a I 



Broos et al. ( 2007 1 in M 17, and similar 



in NGC6357, 
studies. 



5.2. Galactic & Extra- Galactic Gontaminants 

The population of X-ray sources detected in ACIS im- 
ages of Galactic massive SFRs is a superposition of the 
young cluster members, field stars in the Galactic Plane, 
and extragalactic background. In order to examine how 
many X-ray sources toward RCW49 are intrinsic to the 
cluster, we compare the RCW49 popul ation with that 
in a control field of the Galactic Plane (lEbisawa et al. 



[20051). 



he control field, the "Scutum region" at {I, b) ^ (28?5, 
0?0), is an ideal place for our comparison study. First, 
the region was covered by X-ray and NIR imaging stud- 
ies with a sensitivity similar to RCW49. Second, it does 



not harbor a young stellar cluster and the X-ray popula- 
tion there is a typical Galactic Plane population super- 
posed on obscured extragalactic contaminants. Third, 
the extinction in the two regions is similar. The typical 
A^H in RCW49 is derived to be W^-W^-^ cm^^ from 
Figure [sj whic h is consistent with the R OSAT (Belloni 



Sz Mereghetti 19941 and NIR extinction (Ascenso et al 



2007| measurements. The Scutum region shows a bi- 
modal N-ii distribution unlike RCW49 (Fig.js]), but their 
median values are the same (10^^'^ cm~^). in any case, 
the extinction correction is not very important in the 
hard X-ray band, where the opacity is less than one. 

We use the radial profile of X-ray surface number den- 
sity (SND) in RCW49 to argue that a large fraction of 
our X-ray sources consists of the intrinsic RCW 49 popu- 
lation. First, the radial SND profile of all X-ray sources 
(Fig. |6^) shows a centrally peaked shape. The peak at 
R. A.=10:24:01.6 and decl.=-57:45:31, which was derived 
as the mean position of all X-ray sources, is consistent 
with the position of Wd2. Second, the SND of RCW 49 
has an overpopulation compared to the Scutum region 
(Fig. 6j^; dotted histogram) out to ^6'-7'. Within 5' of 
the cluster center where we have SIRIUS coverage, the 
number of X-ray sources is 374 while the projected num- 
ber of the Scutum region is ^69, indicating that ^82% 
of the X-ray sources are intrinsic to RCW 49. 

Our claim that the X-ray population suffers little con- 
tamination by non-cluster-members is more evident for 
those with stellar (optical, NIR, and MIR) identifications 
(Fig. l6b). We have 330 identified X-ray sources within 
5' of the Wd2 center, while the projected number in the 
Scutum region is ~47, which accounts only for ~14% of 
the RCW 49 count. 

This overpopulation holds in the entire X-ray flux 
range of the Ghandra study. In the comparison of the 
hard-band log N — log S relations (Fig. [71, we find the 
total and identified X-ray populations of the RCW 49 re- 
gion (thick histograms) are well above those of the Scu- 
tum region {broken histograms). 

5.3. Distance Constraint from X-ray Data 

The distance to RCW 49 is controversial despite in- 
tensive studies across many wavelengths (§ [2|. In the 
present study, we obtained a rich X-ray sample and 
learned that most of the identified X-ray sources are in- 
trinsic to RCW 49. It is therefore worthwhile to attempt 
to constrain the distance with th e new X-ray da t a set. 

We use a method proposed in Getman et al. (2007b I 
to constrain the distance based on X-ray data. In this 
method, a constraint comes from the dependence of 
X-ray luminosity on the mass of late-type pre-main- 
sequence sources. The relation between these two stel- 
lar quantities is established consist ently in several youn g 
star clusters, including Ceph eus B (|Getman et al.(2006| , 
the Orion Nebula Cluste r ( Getman et al. ' 200d1 i'a- 



ble 6; P7eibisch et al.|2005 1, the I 'aurus Molecular Oloud 
jCiidel et al.p007[ ), and CG12 ( |Getman et aPpOOTb l. 

Note that the X-ray luminosity versus mass relation is 
not degenerate with respect to distance because pre- 
main-sequence X-ray luminosity and mass depend dif- 
ferently on the assumed distance. 

First, we estimate the stellar mass of X-ray sources, 
for which NIR photometry is available (Table |3|. We 
construct a J/(J~H) color-magnitude diagram assum- 
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Fig. 4. — Close-up views of Wd2. (a) SIRIUS Ks-ba nd image. Th e positions of ACIS (crosses), SIRIUS (pluses), IRAC (diamonds), 
and 2MASS (squares) sources are shown. O-type stars in |Rauw et al. (2007) are also shown by thick squares with their names (the prefix 
"MSP" is omitted), (b) Chandra image. The color coding is the same as in Figure [T] (a). The positions of Chandra sources are shown by 
crosses with their sequential numbers (Table [ij . 
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Fig. 5. — Histograms of the A^h values of X-ray sources derived 
by spectral fits (Table [2]l. The solid and dashed histograms are for 
the RCW 49 and the Scutum regions, respectively. 



ing various distances from 0.5 to 8 kpc and derive the 
mass by comparing the d ereddened J-band magnitudes 




to the Siess et al. (2000) model. Then, we extract the 
X-ray sources with an estimated mass of 2.0-2.7 Mq 
and derive the median value of their logarithmic hard- 
band absorption-corrected X-ray luminosity (log L\^ c in 
Table |2| in order to minimize the extinction effect. The 
mass range is limited by the degeneracy in the mass esti- 
mate above ^^2.7 Mq and by our X-ray sensitivity below 
~2.0 Mq. 

Figure [8] shows the median X-ray luminosity at each 
assumed distance, in which the estimated mean X-ray 
luminosity increases monotonically as the distance in- 
creases. We compare the median X-ray luminosi ty in 
the same mass range of the Orion Nebula Cluster (Get- 
man e t al. 2006), for which the distance is reliably mea- 
sured. We found that the distance is constrain ed to be 
2-5 kpc, which is consistent with the claim by (Ascenso 



2 4 6 8 

Radius (arcmin) 

Fig. 6. — Radial SND profile of (a) all X-ray sources, (b) sources 
with optical, NIR, or iVIIR counterparts, and (c) those without 
counterparts. Solid and dashed histograms indicate the data from 
RCW 49 and the Scutum region. The center of the profiles is the 
mean position of all X-ray sources in RCW 49 and the same detec- 
tor position for the Scutum region. The uncertainty is estimated 
by V~N, where N is the number of sources in each bin. 



et al.||2007[ ) but not with ( |Rauw et alT]|2007[ ) 



5.4. Nature of X-ray Cluster Population 

Sources of different natures have different X-ray and 
NIR brightness and colors. In massive SFRs, where the 
mass range of stellar constituents is wider than low-mass 
SFRs, the brightness is of prime importance for estimat- 
ing the nature. We employ X-ray versus NIR brightness 
plots (Figs. [9k and l9b) to estimate the nature of X-ray 
cluster menibers of RCW 49, aided by the conventional 
NIR color-magnitude diagram (Fig. 11 1). 

Figure [9] (a) shows the observed Ag-band magnitude 
{Fks) versus hard-band X-ray flux (F^) of all the X-ray 
sources with Kg-hand identification. Figure [9] (b) is re- 
stricted to a sample with X-ray spectral fits, snowing the 
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Fig. 7.— logAf - logS plots of ROW 49 (solid) and the Scutum 
region (broken). Hard-band absorption-uncorrected flux was used. 
Histograms are constructed for all sources (black) and those with 
stellar counterparts (red), respectively. 
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Fig. 8. — Median X-ray luminosity versus assumed distance 
of RCW49 sources. The logarithmic median value of hard-band 
absorption-corrected X-ray luminosities in the mass ranges of 2.0- 
2.7 Mq is plotted at each assumed distance. The error bars show 
the median absolute dispersion (MAD) divided by the number of 
samples (7-46 sources) at each bin, which represent the uncer- 
tainty of median Lx determination. The dotted and dashed lines 
indicates the median Lx and its MAD divided by the number of 
samples fro m the Orion Nebula Cluster in the same mass range, 
respectively jGetman et al.|2006[ l. 

extinction-corrected absolute i^g-band luminosity (Lks) 
versus hard-band X-ray luminosity {L\^ c)- The constant 
flux ratio lines between the Kg-hand and the hard X-ray 
band {dotted- and- dashed lines in Figure [9^ and b) are a 
proxy for the canonical X-ray versus bolometric luminos- 
ity, which is independent from our distance estimate. 

It is easily recognized that the sources form several 
groups in both plots. The most populous group cen- 
tered at (i^h, Fks) (10-l4■^ 10-12-7) ergs s"! cm'^ in 
Figure [9] (a) and (Lh.c, Lks) ~ (lO^^, lO^S-S) ergs s"! 
in Figure ^ (b) are 1.0-3.0 Mq low-mass pre-main- 



sequence sources in RCW49 for the following reasons. 
Their iTg-band brightness (13-16 mag) is consistent with 
the NIR brightness of pre-main-sequence sources in this 
mass range. Their hard-band X-ray luminosities 10^*^- 
j^qSi.s gj,gg g-i g^-^^ flaring are typical of this class o f stars 
in ne ar-by regions (e.g., ,Wolk et al . 2006^ Getman et al. 
2006| ). The ratio of X-ray to JN IK "brightness of 10"^- 
10~ is consistent with the fractional X-ray to bolometric 
luminosity of 10~''-10~^ common among these sources. 
The _Ks-band bolometric corrections for AO V, FOV, and 
GOV stars are -3.75, ~2-85, a nd -2.25 mag ( |Tokunaga 
2000) iDriUing fc Landolt||2000| . 

!^ number of sources deviate in the upward and 
rightward directions of the primary group in Figure [9] 
Among them, two Wolf-Rayet stars occupy the top 
right end of all sources with {Fu, F^s) ~ (lO"^^-*, 
2^0-10.2) gj,gg g-i Figure [9^ (a) and (ih.c, Lks) 

~ (10^3, 10^6) ergs s-i in Figure [9f(b) . 

The sources in between the late-type pre-main- 
sequence and Wolf-Rayet star groups, both in X-ray 
and NIR brightness, are early-type main-sequence stars. 
They have a large dispersion in the horizontal direction 
in Figure |9] indicating that sources with a similar NIR 
brightness can differ significantly in hard X-ray bright- 
ness. We first extract early-type star candidates based 
on the NIR color-magnitude diagram (Fig. ITok). X- 
ray sources with dereddened iT-band magnitude brighter 
than a B2V star are considered to be early-type can- 
didates and are labeled "ET" in Table [3l A total of 
43 candidate stars include 13 spectroscopically-identified 
0-type stars from preceding studies ( ,Moffat et al. 1991"] 
Rauw et al.|2007 ). These early- type candidates are plot- 
ted by ditt'erent symbols based on their X-ray spectral 
hardness. 

The upward and downward triangle sources, which re- 
spectively indicate hard and soft X-ray emission by X-ray 
spectroscopy, occupy different regions naturally reflect- 
ing the fact that harder O stars are brighter in the hard 
X-ray band. The two groups appear to form two separate 
clustering in Figure [9] and have different X-ray to NIR 
brightness ratio. These suggests that they have different 
X-ray production mechanisms. Eight spectroscopically- 
identifled sources (HD 90273, MSP 151, 171, 182, 183, 
199, 203, and 263) are in the soft 0-star group with 
Lh,c/LKs ~ 10"'', while three (MSP 18, 167, and 188) are 
in the hard 0-star group with L^^c/Lks = lO^'^ — 10^^. 
Source number 29 (ACIS 29) may belong either to the 
hard 0-star or Wolf-Rayet star group, but the nature of 
this source is unknown without any spectroscopic stud- 
ies. We note that the source has infrared excess emission 
(Table [3| and is an off-cluster source at ~4f5 away from 
the center, which is similar to WR20a and MSP 18. 

The soft O stars share common characteristics with 



tradit ional X-ray-emitting O stars (Berghoefer et al. 
19971 for having soft and constant X-ra y emission due 
to mt ernal shocks in their stellar winds (Lucy k. White 
1980|. The hard O stars are among the growin g clas s 



of new-type O stars with hot thermal emission (§ |5.7.2 1 
The three hard O stars (MSP 18, 167, and 188) areTrlght 
enough for detailed X-ray spectroscopy, and we discuss 
them in 



5.7.2 



5.5. Spatial Distribution of X-ray Cluster Members 
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Fig. 9. — Scatter plots of X-ray and NIR brightness of NIR-identified X-ray sources, (a) Hard-band X-ray flux (Fh) versus Ks-band flux 
(Fxs) or magnitude (m/s-s) for 339 X-ray sources with iCg-band counterparts, (b) Extinction-corrected hard-band X-ray luminosity (Lij (,) 
versus i^s-band luminosity (Lks) or absolute magnitude (M/^s) for 124 X-ray sources with Nn estimates in Table[2] Sources are NIR-loud 
toward the top left and X-ray-loud toward the bottom right of the figures. Wolf-Rayet stars are shown by large diamonds. Sources brighter 
than B2V in NIR are shown by upward triangles for bright (more than 50 counts) and hard (best-fit temperature larger than 2 keV), 
downward triangles for bright and soft, and squares for faint (less than 50 counts) sources. Sources fainter than B2 V in NIR are indicated 
by circles. Those with X-ray variabi lity ar e marked with filled symbols. Spectroscopically-identified early-type stars are labeled. The prefix 
"MSP" is omitted for Moff at et al.| ( |1991| l sources. ACIS 29 indicates the source number 29 in Tables Q[2j andjs] The dashed-and-dotted 
lines indicate the iso-x/y lines from 10~' to lO'^ (the logarithmic values are given beside each line). 
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Fig. 10. — NIR diagrams of RCW 49. SIRIUS sources with X-ray and MIR excess emission are respectively shown by blue filled and 
red open circles. X-ray sources with 2MASS but without SIRIUS photometry are plotted by blue filled triangles. 2MASS and SIRIUS 
colors are transformed into th e California Insti tute of Technology color sy stem to compare with the theoretical calculations using the 
formulae given respectively by ICarpe nter] l (2001| andjNakajima et al.' ('2005jl. (a) Color- magnitude diagram. A 2 Myr isochrone curve of 
pre-main-sequence stars (thick solid curve J is from [Barattc et al. (1998 M < 1.4 Mq) and Siess et al. (2000,;. 3.0 Mq > M > 1.4 Mq). 
The reddening lines of 1, 2, and 3 Alp, sources are snown with solid lines. A color-magnitude curve of main sequence stars (thick dashed 
curve) is from Tokunaga (2000); Drilling & Landolt (2000). The reddening lines of 09V and B2V stars are shown with dashed lines. 
The SIRIUS A -band magnitude lim it is in dicated by the dotted line. Spectroscopically-identified early-type stars are labeled. The prefix 
"MS P" is omitted for [Moffat et aT] (|1991[ l sources, (b) Color-color diagr am. The intrinsic colors of dw arfs and giants (thick curves) are 
from |Tokunaga| ( |2000[ ) , the classical 1 laur i star locus (thick line) is from [Meyer, Calvet, &c Hillenbrand| ( [T997j >, and the intrinsic colors of 
Herbig Ae Be stars (dotted curve) are from |Lada &: Adams| ( |T992[ |. 
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We can individually identify cluster members based on 
their X-ray identification with a low level of contami- 
nation. Unlike optical and infrared imaging, the X-ray 
detections are not affected by non-uniform and intense 
diffuse emission. These two features of X-ray imaging 
studies are particularly important outside of the cluster 
core, where the optical and NIR SNDs of cluster members 
fall to a comparable level to those of irrelevant sources. 
Using the X-ray SND profiles, we obtain a sensible esti- 
mate of the spatial scale and distribution of cluster mem- 
bers away from the cluster core. 

In the X-ray SND profile (Fig. [6^), we see two inter- 
esting features. One is that the profile shows an overall 
decline, reaching the background level at 6'-7'. The X- 
ray declining profile is smoother and more extended than 



the N IR profile constructed from 2MASS (Ascenso et al. 
2007 Fig. 6). We obtained a larger estimate ot the clus^ 



ter size b y twice than the NIR estimate of 4^1 (Ascenso et 
ar]|2007 ). Our estimated radius is ~8 pc at the assumed 



Our estimated radius is 
distance, which is far larger than the Orion Nebula Clus 



ter an d somewhat larger than NGC6357 (Wang et al. 
2007al) 



i'he other feature is a deviation from the monotonic 
trend at ~4'. A similar excess can also be found in the 
identified and unidentified X-ray source profile (Figs, pp 
and c). A ^4' radius partial ring around Wd2 is the 
most conspicuous f eature both in the radio and MIR im- 
ages (Fig . [Ill; iWhiteoak fc Uchida|[l997l p^rchwell et| 
al. 20041, where gas and dust are compressed by winds 
and radiation from massive sources in the central OB as- 
sociation. The possible peak at ^4' in the X-ray profiles, 
although the significance is low, suggests that a second 
generation of triggered star formation is present in this 
ring. A deeper X-ray image would confirm this possibil- 
ity and identify members of the triggered population. 

We identified thirty new OB star candidates based on 
the NIR photometry of X-ray sources (§ |5.4[ ). It is inter- 
esting to note that many of them are located outside of 
the cluster core. In addition to the three off-core early- 
type stars (WR20a, 20b, and MSP 18 at -0^6, ~3f7, 
^I'.O from the Wd2 center, respectively; Fig. Tk) consid- 
ered to be physically associated with RCW49, we have 
14 early- type candidates beyond 3' of the cluster center. 
One of them (HP 90273) is a spectr oscopically-identified 
O star (Cannon & Pickering 1919). This implies either 
that these are run- away O stars as suggested for MSP 18 
( Uzpen et al.||200^ l or that the cluster was formed with- 
out primordial mass segregation. We note that the OB 
stars in the Rosette Ne bula Cluster similarl y do not ex- 
hibit mass segregation (Wang et al. 2007b). Follow-up 
studies are mandatory to conhrm the early-type nature 
of these off-core candidates and to understand the cause 
of their large spatial distribution. 

5.6. Discrimination of Intrinsic Cluster Members 

We detected 10,540 and 9,768 sources in the SIRIUS 
images on and off of the RCW49 region (§ |3] hereafter 
called the object and the control fields). This indicates 
that the NIR sample is seriously contaminated by field 
stars. Individual cluster members, most of which are low- 
mass pre-main-sequence sources, are traditionally iden- 
tified by infrared excess emission from their circumstel- 
lar disks and envelopes. MIR exc ess studies are mor e 
sensitive than NIR excess studies (Haisch et al. 2001). 



Whitney et al. (2004) identified 256 cluster members in 
the Chandra field from the MIR excess signature based 
on the MIR color-color diagrams using Spitzer /IRAG. 
However, the infrared excess is only sensitive to classi- 
cal T Tauri stars and protostars with rich circumstellar 
matter. Weak-line T Tauri stars are usually missed. On 
the contrary, X-ray emission in the >1 keV band does 
not significantly depend on whether pre-main-sequence 
sources have disks or not (Feigelson et al. 2007). We now 
use the X-ray emission as a signature of cluster member- 
ship and examine how the Chandra data complement the 
Spitzer data. 

For this purpose, we construct three ifg-band lumi- 
nosity functions (KLFs). One is the intrinsic KLF of 
RCW 49 obtained by subtracting the KLF of the control 
field from that of the object field. The other two are 
the KLFs of SIRIUS sources with X-ray and MIR excess 
emission signatures. We examine what fraction of the 
first KLF is recovered by each and the combination of 
the latter two KLFs. 

Figure [TT] (a) shows the KLFs in the object and the 
control fields {solid and dashed histograms), while the 
difference of the two is plotted by a solid histogram in 
Figure 11 (b). Sources in the two fields have slightly 
different characteristic extinction, for which we compen- 
sate by shifting the object KLF brightward by 0.1 mag 
based on the following argument. The difference between 
the mean NIR colors in the object and control fields 
represents the difference in the extinction. In order to 
avoid the effect of excess emission seen among pre-main- 
sequence sources, we use only the {J-H) color of sources 
without any evidence of infrared excess. The mean (/- 
H) values in the object and control fields are 0.94 and 
0.83 mag with ^0.01 mag uncertainty. The difference 
can be converted to AAk 
Benjamin et al.||2003 l 



0.1 mag (Tokunaga 2000 



l^'irst, we compare the intrinsic KLF of RCW 49 with 
each of the KLFs of sources with X-ray and MIR excess 
emission, which are respectively shown with triangles and 
diamonds in Figure [lT| (b) . Each of the X-ray and MIR- 
excess KLFs represent part of the intrinsic KLF. The 
recovery rate drops at > 14 mag because the Chan- 
dra and Spitzer sensitivities do not match that of SIR- 
IUS. Also, the Spitzer image suffers confusion in the most 
crowded part of the cluster. The recovery of cluster mem- 
bers by the X-ray KLF is also expected to drop in the 
magnitude range for late B and A type sources, because 
some of them are intrins ically X-ray faint unlik e earlier 



and later type sources (Preibisch et al 
sources would appear at 



2005D . These 

12 mag at the assumed dis- 



tance. The smaller X-ray recovery rate at 12.0-12.5 mag 
may be attributable to this effect (Fig.[TT|D). 

We nonetheless do not attempt to correct statistically 
for these types of incompleteness because we are inter- 
ested in how many of the intrinsic cluster members are 
individually identified by the current data set. The total 
counts of sources in the intrinsic KLF at K<14.0 and 
<14.5 mag are 221 and 358. Respectively, 125 (57%) 
and 179 (50%) of them are recovered in the X-ray KLF, 
while 84 (38%) and 98 (27%) are in the MIR excess KLF. 
Each of the X-ray and MIR excess KLFs represent about 
30-60% of the intrinsic RCW 49 KLF. 

Second, we compare the intrinsic KLF with the com- 
bination of X-ray and MIR excess KLFs. The KLFs of 
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Fig. 11. — (a) KLFs of the object and control fields from SIR- 
lUS (solid and dashed histograms). The KLF of the object field 
is shifted leftward by 0.1 mag to compensate for extra reddening 
relative to the control field, (b) Intrinsic KLF of RCW 49 (solid 
histogram) derived as the subtraction of the KLFs of the object 
and control fields. The Poisson uncertainty is given by error bars. 
The KLFs of SIRIUS sources with X-ray and MIR excess (ME) 
identifications are given by triangles and diamonds, respectively. 
The KLFs of SIRIUS sources with either of and both of the two 
identifications (the union and intersection) are given by dotted and 
dashed histograms, respectively. 

the union and the intersection sets of X-ray and MIR ex- 
cess sources are shown by dotted and dashed histograms 
in Figure 11 (b). The union set recovers 90% (200) and 



74% (266jof the intrinsic cluster members at if <14.0 
and <14.5 mag, respectively. The intersection set is 
very small. These indicate that X-ray and MIR excess 
emission selections operate in a complementary fashion 
in identifying cluster members and that the combina- 
tion of the two recovers a large fraction of the mem- 
bers. The complementarity probably comes from the 
fact that richer circumstellar matter around pre-main- 
sequence sources increases the MIR excess detection rate 
due to stronger emission and decreases the X-ray detec- 
tion rate due to larger extinction. This i s be st examined 



in the NIR color-color diagram in Figure 10 (b). SIRIUS 



sources with X-ray and MIR excess signatures (blue and 
red dots) are segregated in the plot. The X-ray emission 
preferentially detects sources with no or smaller NIR ex- 
cess. 

5.7. X-ray Emission from Massive Stars 

We present the results of detailed X-ray spectroscopy 
for five exceptionally bright X-ray sources with more 
than 200 counts. Two of them are Wolf-Rayet stars 
(WR20a and 20b) and the rest are 0-type stars (MSP 18, 
167, and 188). All these sources are located at the top 
right end of X-ray and NIR brightness plots (Fig.|9]). 

5.7.1. Wolf-Rayet stars 

The X-ray spectra of WR 20a and WR 20b are shown 
in Figure 12 (a) and (b), respectively. Both of them are 
characterized by emission lines from highly ionized ions 
such as Mg, Si, S, Ar, and Ca, and a hard tail up to 
~5 keV. The Si XIII Ka emission as well as the hard tail 
indicate that the emission is more likely from a plasma of 



multiple temperatures than from that of a uniform tem- 
perature. We present the results of a two-temperature 
plasma fit in Table |4] In the fitting, the abundances of 
conspicuous lines (Si and S) were thawed to obtain the 
best-fit values. The abundances of other elements were 
fixed to be 1 solar. It is known tha t hydrogen is heavily 



depleted in Wolf-Rayet stars (e.g., van der Hucht et al. 
1986 Morris et al. 2000), but we used the cosmic abun- 
dance ( p\.nders G rcvcssc 1989) as we do not know the 
level of hydrogen depletion for these stars. Once it is 
known, the abundance values presented here should be 
decreased accordingly. 
WR 20a is an eclipsing binary of two Wolf-Rayet stars 



of spectral type WN6ha (Bonanos et al 



2004 



Rauw et 



al.''2005). Using the best-fit orbital solution by Bonanos 
et al. (2004), the Chandra observation covers the orbital 
phase from 0.36 to 0.48. Because of the uncertainty in 
the orbital period determination of 3.686 ± 0.01 days, 
the Chandra phase range has uncertainty from 0.18-0.30 
to 0.55-0.67. We compare the X-ray light curves and 



the eclipsing optical light curve (Fig. 13) with the aim to 



obtain a hint in the X-ray production mechanism of this 
source. One artificial effect is introduced by the fact that 
the source was observed in a gap between ACIS 12 and 13. 
The source oscillated across the gap at a period of 1000 s 
due to the dithering motion of the telescope. About 60% 
of the photons were lost in the gap. We binned the light 
curves with a binning size of 1000 s to cancel this effect. 

No apparent variation wa s fo und in either the aver- 
age X-ra y count rate (Fig. 13 j) or the hardness ratio 
(Fig. 13"). In contrast, a clear decline is seen in the I 



band magnitude (Fig. 



13 i) during the phase of about 
0.3-0.7. The Chandra phase coverage is uncertain from 
0.18-0.30 to 0.55-0.67. When its lower limit is adopted, 
the Chandra observation was conducted out of eclipses. 
When other values are adopted, it was conducted during 
an optical eclipse. In this case, the lack of the eclipse in 
the X-ray light curves indicates that the X-ray emitting 
region is far enough from the stellar surface so as not 
affected by the photosph eric eclipsing. From the orbital 
solution by Bonanos et al. (2004), the radii of the two 
^19Rq and asin(i) is ~53i?0, where a is the 
i =74.5° is the inclination 
The system is 
At 



stars are 

semi-major axis of the orbit, 
angle, and Rq is the radius of the Sun. 
observed close to the edge-on in the line of sight, 
the middle of the eclipse, the two stars are separated by 
2a cos (z) ~ 29i?0. Therefore, the X-ray emitting region 
is extended to at least a few stellar radii. 

5.7.2. 0-type Stars 

At least two out of three bright O stars show emission 
line features of highly ionized Mg, Si and S, indicating 
their thermal plasma origin. We fit all three spectra with 
a thin-thermal plasma model with variable abundance of 
elements with conspicuous features. Table |4] gives the 
results. All of them are characterized by a hard spectrum 
with the best- fit plasma temperature of >1.5 keV and 
luminous X-ray brightness of 2-8x10'^^ ergs s~^. 

Such hard and luminous X-ray-emitting O s tars have 
been report ed in some other massi ve SFRs ("Albacetel 
Colombo et al. 2003: Rho et al. 2004: S kinner et al. 200 ~ 
Broos et al. 2007: Wolk ct al. 2006). The X-ray proper 
ties of these sources arc quite ditt'crcnt from the classical 
view of main-sequence 0-type stars characterized by soft 
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Fig. 12. — X-ray spectra of bright sources. The source identifications and the models employed to fit the spectra are given in Table |4] 
Grouped data with errors are shown in the upper panels over-plotted with the best-fit model convolved with mirror and detector responses 
(solid line). The energy of the Ka lines of some representative highly ionized ions are shown. The lower panels show the residuals of the 
fit. 

ries" if this interpretation is correct. Yet another idea 
is that inverse Compton emission of UV phot ons from 
these O stars accoun ts for their hard spectra (Albacete 
Colombo et al.|[2003 ). The X-ray production mechanism 
of such curious O stars is not yet certain, but enriching 
the sample using plots like Figure |9] may lead to a better 
understanding of the mechanism. 

6. SUMMARY 

We observed the Galactic massive star-forming region 
RCW 49 and its central OB association Wester lund 2 us- 
ing Chandra/ AGIS, IRSF/SIRIUS, and Spitzer /m.AG. 
The main results are summarized as follows. 

1. 468 X-ray sources were detected in a ~40 ks ACIS 
exposure covering a ~17'x 17' field. Previous X-ray 
studies detected less than 10 sources in the studied 
field and almost all the X-ray sources reported here 
are new detections. The results of X-ray photo- 



thermal spectra of <1 ke V and the X-ray to bolo metric 
luminosity ratio of ~10~^ ( Berghoefer et al.|l997 l. Hard 
X-ray-emitting O stars in massive SJ^'Ks should be a new 
distinct class of X-ray sources with a different X-ray pro- 
duction mechanism. 

Several interpretations have been proposed for such 
hard X-ray emission fr om O stars. One is magnetically- 
confined wind sh ocks ( Babel & Montmerle 1997: Schulz 
[eraH2003| |Gagn^ et aF ^OOS; Stelzer et al. 2005), where 
magnetized winds produce hard X-ray emission in large- 
scale shocks near the stellar surface. Another is similar to 
the X-ray emission in WR binaries; i. e., winds from each 
component of a close binary of two 0-type stars collide 
with each other to make h ard X-rays dZhekov & Skinner 
2000[ [Skinner et"al]|2007[ ). The hard^X-rays might be 
from the mag netic reconnections between the two binary 
components ( Schulz et al. 2006 1 . The binary nature of 
the three 0-type stars in ROW 49 is not confirmed yet, 
but these sources might be "X-ray-spectroscopic bina- 
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Fig. 13. — Light curv es of the binary system WR 20a. (a) /-band 
magnitude of one orbit | |Bonanos et al.|2004[ |. (b) X-ray count rate 
in the 0.5-8.0 keV band and (cj the hardness ratio (HR) in the 
phase 0.36 to 0.48 (the range is indicated by two dotted Hues in a; 
the un certainty is about ±0.2). We used the ephemeris by |Bonanos| 
|et al. 1 120 04) with an orbital period of 3.686 days and a JuTian date 
ongm o! 2,453,124.569 day HR here is defined as {H-S)/(H-|-S) 
where H and S are the count rates in the hard (2.0-8.0 keV) and 
the soft (0.5-2.0 keV) bands, respectively. 

metric, timing, and spectroscopic analyses are pre- 
sented and summarized in Tables [l] and [21 

2. We obtained SIRIUS images in the J, H, and Kg 
bands of a concentric ~8f3x8'3 field and extracted 
10,540 sources at a 10 a detection limits of ^-^19.0 
(J), -17.8 {H), and ~16.8 {K^) mag. 

3. 379 X-ray sources were identified with optical, NIR, 
and MIR counterparts using NOMAD, 2MASS, 
SIRIUS, and GLIMPSE data. The results are given 
in Table IS 

4. The central OB association Westerlund 2 was 
resolved in the X-ray image for the first time. 
X-ray emission was detected from all of the 
spectroscopically-identified early-type stars con- 
sisting of two Wolf Rayet stars and a dozen O type 
stars. 

5. About 86% of the X-ray sources with optical or in- 
frared identifications were found to be cluster mem- 
bers by comparing with the X-ray population in 
a control field in the Galactic Plane. The X-ray 
overpopulation against the control field is seen in 
all measured fiux ranges. 

6. We compared the intrinsic if -band luminosity 
function of NIR sources in RCW 49 with the KLFs 
of sources with X-ray or MIR excess identifica- 
tion. We found that 30-60% of the cluster mem- 
bers at Xs<14 mag are individually identified using 
the two observational signatures of cluster member- 
ship. When X-ray and MIR excess signatures are 
combined, about 90% of cluster sources are indi- 
vidually located. The two indicators of the cluster 
membership work in a complementary fashion, in 
which the former excels in detecting weak-line T 



Tauri stars with smaller attenuation while the lat- 
ter predominantly detects classical T Tauri stars or 
protostars with richer circumstellar matter. 

7. A loose constraint on the distance to RCW 49 was 
derived to be 2-5 kpc based on the mean hard-band 
X-ray luminosity of T Tauri stars in the mass range 
of 2.0-2.7 Mq. 

8. The cluster X-ray population consists of low-mass 
pre-main-sequence sources and early-type sources 
based on X-ray and NIR photometry. Late-type 
pre-main-sequence sources, Wolf-Rayet stars, and 
early-type main sequence stars occupy different re- 
gions of X-ray versus NIR brightness plots. Thirty 
new OB candidates were identified based on the 
NIR color-magnitude diagram, several times more 
than the number of optically-identified OB stars. 

9. Using the X-ray surface number density radial pro- 
files, we gave an estimate of the cluster radius to be 
6'-7', or ~8 pc at the assumed distance. A possible 
excess is seen at ~4' of the center, suggesting that 
secondary star formation is taking place at the ring. 
Fourteen OB star candidates are found outside of 
Wd2, suggesting an absence of mass segregation. 

10. Detailed X-ray spectroscopy of the two Wolf-Rayet 
stars was presented. Both sources show two- 
temperature thermal plasmas with characteristic 
temperatures of 0.5-1.0 and 4-6 keV and 1-8 keV 
luminosities of 2x10^^ ergs s~^. No X-ray fiux 
or hardness variation was found from WR20a, an 
eclipsing binary, suggesting that the X-ray emission 
arises far from stellar surfaces. 

11. Detailed X-ray spectroscopy of three luminous O- 
type stars were presented. They are characterized 
by thermal spectra with a temperature of >1.5 keV. 
These are among the growing samples of hard X- 
ray emitting O stars recently found in some massive 
star-forming regions. Early-type stars appears to 
have two sub groups, hard and soft, based on the 
X-ray spectroscopy. 
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